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The problem of a spatially correlated noise affecting a complex system is studied in this paper. We 
present a comprehensive analysis of a 2D model polymer chain, driven by the spatially correlated 
Gaussian noise, for which we have varied the amplitude and the correlation length. The chain 
model is based on a bead-spring approach, enriched with a global Lennard- Jones potential and 
angular interactions. We show that spatial correlations in the noise inhibit the chain geometry 
dynamics, enhancing the preservation of the polymer shape. This is supported by the analysis of 
correlation functions of both the module length and angles between neighboring modules, which 
have been measured for the noise amplitude ranging over 3 orders of magnitude. Moreover, we 
have observed the correlation length dependent beads motion synchronization, and the spontaneous 
polymer unfolding, resulting from an interplay between chain potentials and the spatially structured 
noise. 
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I. INTRODUCTION 

The understanding of diffusion in complex media is 
crucial for both the modeling of conformation transitions 
in biomolecules and the intracellular transport. It is also 
well known that various systems organize spontaneously 
in a response to random forcing pQ and that the introduc- 
tion of temporal correlations into the noise can lead to 
synchronization effects [2 j. A well-established framework 
to simulate these phenomena is provided by Langevin 
equations, which introduce the concept of stochastic force 
mimicking the molecular collisions [3]. An important 
advance in this formalism has been the introduction of 
Generalized Langevin Equation (GLE), which reproduce 
the anomalous diffusion thanks to the time-correlated 
stochastic force and the corresponding integral memory 
kernel, which represents the friction [4]. Recently, S. C. 
Kou in [4] has derived the GLE from a microscopic model 
of a particle coupled to a large number of oscillators, 
thus showing that the particle-environment interaction 
is essential for the occurrence of temporal correlations in 
thermal noise. However, it is remarkable that this the- 
ory explains solely the temporal aspect of diffusion, while 
little work has been done to understand its spatial coun- 
terpart. This has led us to investigate the problem of a 
spatially correlated noise affecting a complex system. 

The collective media behavior which is random but 
characterized by a certain correlation length A occurs at 
the length-scale of micrometers in the context of hydro- 
dynamic interactions, e.g. in colloid sedimentation [3 [6] 
or in the study of the active particle motion [7 j. However, 
the spatial correlations at the lower length-scale play a 
fundamental role in the theory of phase transitions [8] 
among which the liquid-glass transition is of special in- 
terest. During this transition, the particles suffer a dra- 
matic drop of mobility without the emergence of struc- 
tural ordering [9 . This phenomenon has been intensively 
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researched on for past two decades, and, according to ex- 
tensive simulations [10] [11], it is characterized by the 
occurrence of the spatial correlations in particles motion 
[TT] , which is recognized as the formation of the different- 
sized clusters p~0j[T2]. Choosing a single moment in time, 
one could interpret these clusters as a source of a distur- 
bance which is random at the large length scale (^> A), 
but ordered at the length scale of A. Figure [I] illustrates 
this idea. The temporal evolution of this system is still 
indeterministic, as it 'randomly reorders'. We propose 
that this behavior could be imitated by the spatially cor- 
related noise, which is affecting a subsystem, in our case, 
a model 2D polymeric chain. 

We have simulated the chain based on the bead-spring 
approach under the forcing of Spatially Correlated Gaus- 
sian Noise (SCGN) for which we have varied the corre- 
lation length and the amplitude. Our previous findings 
regarding the stiffening of the chain under the SCGN, 
shown with the aid of the reduced dynamics, have been 
published in [15] . However, our further investigation into 
this system, which involves the extension of the parame- 
ters' range and the measurements of chain characteristics, 
has revealed several new effects, namely: beads motion 
synchronization, increased time correlation of both mod- 
ule length and angles between modules, the inhibition 
of the average module length growth and, most notably, 
the chain unfolding induced by the increased correlation 
length. 

Our simulations are related to the actual physical sit- 
uation by the choice of A. Unfortunately, currently there 
are few experimentally accessible quantities that describe 
the collective molecular behavior in the vicinity of glass 
transition and can be measured for the variety of tem- 
peratures P3] . One of these parameters is the number of 
cooperatively rearranging molecules [15] , which has been 
reported to rise from 1 in liquid phase to approximately 
10 in the glass phase [15]. Additionally, these results 
are qualitatively similar for the different chemical com- 
pounds [15]. On the other hand, the direct measurements 
of the correlation length are scarce and limited to spe- 
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cific experimental setup, as e.g. in [16 which reports A to 
be of the order of 2-4 molecule's diameters. These mea- 
surements suggest that A covering up to 5 chain nodes is 
physically meaningful. 

The paper has following structure: in section [IT] the 
methods of the SCGN generation are introduced, in sec- 
tion [ITT] we propose the equations of motion and the cor- 



relation function, in section IV we present our polymer 
model, sectio n |V| briefly discusses simulation methods, 
sections from | VI| to [IX| present the results regarding each 
effect with an interpretation, in section [X] we summarize 
our findings. 



II. MULTIPLE CORRELATED GAUSSIAN 
VARIABLES 

The generation of multiple correlated Gaussian vari- 
ables is a central problem in the simulation of SCGN 
driven systems, therefore we shall outline here the basic 
algorithm. 

Lets assume that we have two real vectors of random, 
zero-mean Gaussian variables, namely = (£i, . . . , £/v) 
and ff = (r/i, . . . , t]n), which components satisfy the fol- 
lowing correlation relations: 



(ViVj) 



So 
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Here, Sij denotes the Kronecker delta, and Sij are ele- 
ments of the correlation matrix, defined as: 



(3) 



The matrix S is symmetric and positively definite [17] , so 
it is suitable for Cholesky decomposition [18] , which fac- 
torizes S into a lower triangular matrix L and its trans- 
position: 



S = LL 



(4) 



The vector of correlated variables £ is related to the un- 
corrected vector ff via a linear transformation [19]: 



£ = Lrf 



(5) 



This means that, given a correlation matrix, one can gen- 
erate the correlated Gaussian vector £ simply by sam- 
pling N times the normal distribution to obtain the com- 
ponents of ff and then performing the transformation ([5| . 



III. EQUATIONS OF MOTION AND 
CORRELATION FUNCTION 

Our system is equivalent to an ordered set of N in- 
teracting material points on a plane, enumerated by the 
index i. The position of i-th point (or bead, as we will 



refer to it further) is ff = {x^yij. In order to simu- 
late the trajectory {fiit)}^ of the whole system, we have 
to solve numerically a set of 2N stochastic equations of 
motion: 



mxi + jXi 
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Here, U is the potential energy of the system, which 
we will discuss in detail in the next section. £,(fi) T = 
(£a?(^)j is the two dimensional SCGN, m is a bead 

mass and 7 is a friction constant. In an absence of a 
more relevant theory, we have applied the simplest fric- 
tion model, and chosen 7 to be constant. The differential 
equations (|6| are in principle of the second order, which 
we preserve for generality, but in the course of our sim- 
ulations we have overdamped the system by choosing 7 
to be large enough. 

We assume that the correlation function Sij of stochas- 
tic forces acting on beads i and j should depend only 
on a relative distance between these beads, which is 
fij = \u — Additionally, we also assume that there 
are no cross correlations between x and y components, 
which allows us to reduce the correlation relations to the 
form: 



(€x(fi)€ x (fj)) = (€ y (fi)€y(fj)) = S(rij) 



(7) 



It should be emphasized that the correlation matrix 
S is a dynamical object, and evolves in t as the rela- 
tive distances rij(t) do. The conditions ^ suggest the 
following procedure to integrate equations (JoJ) : once all 
bead's positions {fl(t)}N at some moment t are deter- 
mined, we can calculate the NxN correlation matrix and 
its Cholesky decomposition L; next, according to ([5|, we 
shall use L and two different ff to obtain {^(t^)}at and 
{£ y (fi)}N- Finally, we can use them to perform an in- 
tegration step, which gives {fl(t + At)}w- The repeated 
Cholesky decompositions are the most computationally 
expansive part of our simulations, as the computational 
complexity of this decomposition is 0(7V 3 )[18 . 

Along with conditions (|7|), we assume that the cor- 
relation function Sij is characterized by the correlation 
length A, and it reproduces the standard Brownian dif- 
fusion for A — » [20], so: 
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In the above formula k denotes Boltzmann constant and 
T is a temperature. Taking into account ([7]) and ([8|, 
we chose the exponentially decaying spatial correlation 
function, which resembles the displacement correlation 
function from [11] and [T6] . We also neglect the tem- 
poral correlations, as we are interested in the effects of 
the purely spatially structured noise. Finally, the spatio- 
temporal correlation function reads: 
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Figure 1. A spatially correlated random vectors with corre- 
lation length A, generated on a regular network. For A = 
the pattern is entirely random, but for A = 20 the ordered 
clusters can be noticed. 



a = 2kT denotes the noise amplitude and we will refer 
to it as temperature, as it is proportional to the actual 
physical temperature. 

In order to illustrate how the spatial correlations affect 
the noise pattern, we have applied ([5| and ([9| to generate 
the random vectors on a regular network. A snapshot 
from this simulation is presented in Figure [I] One can 
easily notice clusters of correlated vectors, however, this 
pattern changes dramatically for every new generation. 



IV. THE MODEL OF POLYMER CHAIN 

The polymeric chain is an archetype of many 
biomolecules, thus we have chosen it as a test-object for 
our simulation. Our model is based on the bead-spring 
approach, in which i and i + 1 beads interact with a har- 
monic potential: 

N-l 

2=1 

Every bead is also the source of the Lennard- Jones type 
interaction, which provides excluded volume effect and 
an interaction between the distant tails of the chain: 
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Finally, we introduce a harmonic interaction between 
beads i and i + 2, which resembles angular interactions: 

N-2 

The total potential energy U is equal to: 

U = U R + U^ + U lj 

For e = (no Ulj contribution) the potential energy 
is minimized when beads' positions satisfy: 

- n\ = d 
1^+2 - n\ = l 



(12) 
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Table I. The parameters of the system chosen for simulation. 



In this case, all of minimum energy conformations are 
equienergetic. In fact, unless lo > 2<io, once f\ and r<i 
are chosen to satisfy \?2 — r\\ = do, the 3rd bead can be 
positioned in two ways, so the relation \fs — f\\ = lo is 
also fulfilled. Successively applying the conditions (14) 



to following beads, one can built numerous minimum en- 
ergy geometries. When Ulj ^ the energetic structure 
of the chain become more complex, but if do > (Jlj and 
e c± fc 1? the Lennard- Jones contribution becomes a per- 
turbation. However, the Ulj influence makes the struc- 
tures no longer equienergetic. 

When the chain's energy is not minimized, the dynam- 
ical topography of the potential energy surface depends 
on both potentials' parameters and the local geometry of 
the chain. An effective way to represent snapshots of this 
energy landscape for a single bead is to take into account 
its 4 nearest neighbors. An example of such a landscape 
is reproduced in Figure [2]A. We have chosen the values 
of potential energy parameters (Table |l| such that the 
double minimum structure is distinct and holds for the 
wide range of local conformations. However, this struc- 
ture is extremely sensitive to a single parameter, which 
is the distance lj = |r^+i — Whenever lj > 2do, 

the two minima tend to merge rapidly into a single one, 
positioned in-line with the beads j — 1 and j + 1. This 
is reproduced in Figure [2^. This fact significantly affects 
the high temperature dynamics of the chain, as is show 
in section HXl 



V. SIMULATION 

Applying the classical Runge-Kutta method modified 
for stochastic differential equations [21], we have simu- 
lated the system described by equations d6) with the po- 
tential ( 13 ) and the parameters from Tablellj The number 
of beads has been set to TV = 128, the bead's mass has 
been chosen m = 1 and the friction coefficient 7 has been 
set to 20, which overdamped the system. 

In our research, we have explored three regimes of tem- 
perature. First, we have varied the noise amplitude a 
from to 20 units at the interval of 1 unit, and we have 
increased the correlation length A from to 20 at the 
interval of 5 units. In the second regime, we have in- 
creased a from 25 to 250 at the interval of 25 units, and 
in the third regime we have explored region from 300 to 
1000 units, at the interval of 100 units. For the second 
and the third regime we have varied A from to 50 at the 
interval of 10 units. For each pair of A and a we have per- 
formed 64 runs, starting from different initial positions. 
The initial coordinates has been chosen, so the distance 
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Figure 2. The energetic landscape for a single bead interacting 
with its 4 nearest neighbors. Potential parameters are chosen 
according to Table [I] A) The distance — Tj-\\ — /o, ji 

and j2 enumerate two possible positions of j-th bead that 
minimize the potential energy. U(0, 0) is the height of energy 
barrier and U m in is the depth of the minimum. B) \rj+i — 
fj-i \ = 1.25/o, ji and J2 merge into a single minimum as the 
energy barrier disappears. 



VI. THE BEADS MOTION 
SYNCHRONIZATION 

The introduction of the spatial correlations into the 
noise implies that, at the length scale comparable to the 
correlation length A, the stochastic force vectors have 
similar direction and value. Therefore, one could expect 
that the motion of beads, which relative distance is lower 
than A, will synchronize. This prediction has been fully 
confirmed. 

As the measure of synchronization at a particular mo- 
ment t, we have chosen the normalized product of two 
beads' velocities, distanced by n nodes, which has been 
averaged along the chain: 

N—n _ 

*•<«> - w^) g " (15) 

Here 6i^ +n is an angle between velocity vectors of the z-th 
and z + n-th bead. For each run, we have gathered K n (t), 
which was time-averaged to obtain synchronization factor 
K n . The maximal value of K n = 1 indicates a fully 
synchronized motion, while K n = implies the opposite. 

We have gathered the data for n ranging from 1 to 9. 
A representative sample of our results is show in Figure 
[3] The rise in the synchronization factor K n along with 
increasing A and a = const is evident. Conversely, the 
level of synchronization is almost constant for A = const 
and varying amplitude, which is valid even for tempera- 
tures below cr = 5. For every n, the factor K n grows from 
for A = to the maximal observed value for A = 50, 
which is approximately 0.8 for n = 1 and 0.3 for n = 9. 

A further insight into the synchronization comes from 
the rearrangement of data, so K n is represented as a func- 
tion of n with A and T being parameters. Figure [4] shows 
the qualitative similarity between these data for two ex- 
treme temperatures (a = 1 and a = 1000). To obtain 
a quantitative measure of the decrease in synchroniza- 
tion with the rise in n, we have fitted our data with the 
exponential decay model: 

K n = A x c~ B - n (16) 



between nearest neighbors has been equal to do? but the 
angle between modules has been chosen randomly from 
tt/2 to 3tt/2. 

The integration step has been set to 1/128 time unit, 
and each simulation lasted 2148 time units. The data 
for the first 100 units has been rejected due to the sys- 
tem thermalization. If not stated otherwise, the data 
has been collected once per time unit. We have gathered 
the data regarding beads synchronization, module length 
and angle between modules. 



This model proved to be an accurate description of data, 
as the coefficient of determination R 2 exceeded 0.99 for 
all fits, except those with A = 0, for which B\=o has no 
physical meaning. 

In Figure |5j we have juxtaposed the values of B\ for 
o~ > 25, at which temperature the behavior of chain is 
noise-dominated. According to this figure, the value of 
B\ is mainly determined by A and decreases when the 
temperature grows by two orders of magnitude. How- 
ever, for A > 30 this fall is rather insignificant, thus, we 
conclude that the noise correlation length is the primary 
factor that influences the effective range of synchroniza- 
tion along the chain. 
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Figure 3. The synchronization factor K n for n — 1,4, 8 as a function of the correlation length A and the noise amplitude a. 
Each column contains the data for the same n in the high temperature regime (upper graph) and the low temperature regime 
(lower graph). 



VII. BEADS MOTION CORRELATION 

The other quantities that are also affected by the pres- 
ence of spatial correlations in noise are: the time cor- 
relation of the modules length and the time correlation 
of the angles between neighboring modules. These two 
characteristics describe the time evolution of the chain 
geometry. 

By a module we understand two neighboring beads, so 
the length of the j-th module, at certain moment £, is 
denned as: 

d j (t) = \f j (t)-f j . 1 (t)\ (17) 

The angle between two neighboring modules is denned 
by the positions of the three following beads: 

iPjit) = L fo_i(t), ^ + i(t)) (18) 

With the beginning at the center of the coordinate sys- 
tem, vectors T{ are equivalent to the coordinates on a 
plane, thus they are applied in the above definition. Ad- 
ditionally, one has to remember that the angle ipj is di- 
rected and varies from 0° to 360° (with 180° indicating 
that three beads are exactly in-line), so the angles have 
to be measured in a unified way along the whole chain, 
conserving the initial numeration of beads. 

The time correlation function of angles has been cal- 



culated in a following way: 

T N-l 

c*{t) = — J2 E mmtk) - <v» 

°* k=0 j=2 

(19) 

Here, we introduce the additional summation over j due 
to the fact, that we have N — 2 angles for a single mo- 
ment t, which allows us to increase statistics and obtain 
a correlation measure for a whole chain, rather than a 
single site. The normalization factor has been chosen 
as: 

T N-l 
k=0 j=2 

which means that C^(0) = 1. Finally, (ip) reads: 

1 T JV-1 

w = tmEE^^ ( 21 ) 

v } k=0 3=2 

In a strict analogy to the angle correlation function 
CV(r), we can define the module length correlation func- 
tion C d (r): 

T N 

C d(r) = ^EE^'^ + T ) " WM^) - < d » ( 22 ) 

A;=0 j—2 

The Cd and (d) are defined similarly to their angle coun- 
terparts. 
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Figure 4. The synchronization factor K n for two extreme 
values of temperature, as a function of the beads' distance n. 
The data has been fitted with the model: K n = A\e~ Bxn . 



An example of collected data is presented in Figure 
[6j In the high temperature regime (approximately for 
a > 200), both C^{r) and Cd{r) are positive functions, 
asymptotically falling from 1 to 0, typical of the stochas- 
tic motion. However, they differ significantly in the low 
temperature regime. While C^(r) preserves its high tem- 
perature profile (but with values much closer to 1), the 
Cd(r) reassembles a linear function, falling below with 
the increase in r . This long-term behavior of the low- 
temperature Cd{r) indicates the domination of the deter- 
ministic motion in this temperature regime. In the con- 
text of the energetic landscape, introduced in section [TV] 
we can suppose that the beads are trapped at the bottom 
of their potential energy wells and perform the damped 
oscillations, slightly perturbed by the noise. Apparently, 
while the beads' motion makes dj oscillate, it barely af- 
fects the angles, so the values of C^(r) are relatively close 
to 1. Additionally, the comparison between C^(r) and 
Cd(r) suggest that the module length behavior evolves 
from deterministic into purely stochastic one as the tem- 
perature grows, while the ipj (t) is of the stochastic nature 
for every a. 

In order to measure the influence of a and the noise 
correlation length A on the C^(t) and Cd(r), we have 
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Figure 5. The synchronization decay factor B\ as the function 
of temperature, for a > 25. 



fitted the profile functions with the following model: 

C(t) =ar a +c (23) 

Despite inaccuracy for r — >• and a divergence in the 
low-temperature regime, the power function model pro- 
vides a quantitative information on a and A dependen- 
cies, thanks to the a parameter. The values of ad and 
a^ plotted against the a and A are shown in Figure [7[ 
As expected, for all values of a and A, a is negative, and 
tends to with the decrease in temperature. However, 
for constant A, a^ decreases at a similar pace with the 
growth of (j, while ad varies slowly for the most of the 
temperature range, but jumps rapidly below a = 100. 

The increase in the noise correlation length A affects 
both ad and a^ in a similar way, namely, the larger A, 
the lower \a\ is obtained. This means that the correlation 
functions C^(r) and Cd(r) decrease at slower rate and so 
the dj(t) and ipj(t) vary less rapidly over time. Therefore, 
the dynamics of chain's shape becomes attenuated and a 
current conformation is preserved longer. 



VIII. MODULE LENGTH DISTRIBUTION 

We have also investigated the marginal distribution 
T(d) of the module length d and its temperature evo- 
lution, with and without spatial correlations in noise. 
Taking into account that dj (t) may express an oscillatory 
behavior, we have reduced the time interval between data 
acquisitions to 1/4 of time unit, to avoid synchronization 
effects. The spatial resolution of histograms has been set 
to 0.21 length unit. 

The profile of T(d) proved to be a single peaked dis- 
tribution, concentrated in the vicinity of its mean, with 
slight, but noticeable asymmetry. Therefore, in order to 
describe T(d) we have calculated its first moment, second 
central moment (presented in Figure [8| and the skewness. 
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Figure 6. A representative selection of the module length correlation functions Cd(r) and the correlation functions of the angles 
between neighboring modules C^(r). Profiles has been fitted with function ar a + c. 



-0.5 \ 
-1.0 



" 1000 



var(d) 




Figure 7. The exponents (A) and (B) resulting from 
fitting the power function to correlation profiles Cd(r) and 
C*(t). 



Figure 8. Statistical characteristics of the module length dis- 
tribution r(d) as the functions of the noise correlation length 
A and the temperature a: A) the mean module length (d), B) 
the variance of T(d). 



The skewness grows with <j from to approximately 
0.6 and saturates at this value. Fortunately, asymmetry 
proved to be small enough, so the others parameters are 
still physically meaningful. The dispersion of distribution 
Y(d) grows with the increasing a (Figure [8^3), which is 
expected diffusive behavior, but also the mean distance 
(d) tends to grow (Figure [8^), starting from (d) ~ do. 
This fact, along with the non-zero skewness, indicates 
that the underlying potential is asymmetric, and, indeed, 
the presence of the repulsive Lennard- Jones core provides 
the reflective barrier preventing two neighboring beads 
from closing up. Conversely, the lengthening of dj is still 
possible as the energy well is not so steep for dj > do as 
in the opposite situation. 

The spatial correlations in noise play an inhibitory role 
for the process of the temperature dependent broadening 
of r(d). A ^ slows down the growth of both (d) and the 
dispersion of T(d). This effect can be explained by the 
following reasoning. When thermal bath imposes non- 
correlated, stochastic forces on beads j and j + 1, this 



commonly results in a nonzero relative force stretching 
(or shrinking) the module. However, when the noise is 
spatially correlated, stochastic forces applied to beads be- 
come similar at the length-scale of A, which significantly 
reduces the relative forcing and, in turn, the d is less 
affected by the noise. 



IX. POLYMER UNFOLDING 

The most unexpected effect that stems from the pres- 
ence of the spatial correlations in the noise is the sponta- 
neous linearization of the chain. Having defined the an- 
gles ipj(t) in (18) we have been able to obtain a marginal 
distribution of angles depending on temperature a 

and correlation length A. Similarly to previous section, 
the data has been collected every 1/4 of time unit, with 
the resolution of histogram set to 1°. The representative 
selection covering the entire range of tested parameters 
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Figure 9. The dispersion of distribution as a function 

of the noise correlation length A and the temperature a. 



is presented in Figure [10 



The temperature evolution of distribution gives 
an insight into how the angular degrees of freedom are 
freed with the rise in temperature. Let us analyze the 
A = case, first. For low temperatures (a < 10) we 
obtain a symmetric bimodal distribution, which is in ac- 
cordance with the predictions of the double minimum 
energetic landscape. However, for an extremely low tem- 
peratures (a < 3) one can see 4 distinct peaks, which 
indicates that, probably, there are 2 additional minima. 
We can suppose that they are shallow, as they disappear 
fast with the rise in a. For temperatures from a = 5 to 
a = 13 the increased penetration of the energy barrier 
region is viewable, and the third peak appears exactly 
at ip = 180°. At a = 13 the two peaks indicating en- 
ergy minimums can be no longer distinguished and from 
now on, the shape of distribution gradually transforms 
from a bell-like curve into a triangle profile, which is 
completed approximately for a = 100. Since then, the 
distribution broadens systematically with the increase in 
temperature. 

The introduction of A ^ affects $(V0 in a subtle, but 
remarkable way. In the low noise regime (a < 25) the 
increase in A retards the temperature evolution of ^(VO? 
so the bimodal profile is preserved in a wider interval of 
a. However, when a exceeds 50, the profile transforms 
into a heavy-tailed peaked distribution, much more con- 
centrated in the vicinity of ip = 180° than in the case 
of A = 0. The dispersion of <&(i/)) is a measure of this 
effect, which shows that the higher A, the lower value of 
the $( , 0)'s second central moment. This indicates the 
linearization of chain, and thus, its unfolding. It is illus- 
trated in Figure [9j 

The single-peaked distribution with (ip) = 180° indi- 
rectly suggests that for high temperatures the system se- 
lects the single-minimum energetic configuration similar 
to the one shown in Figure [2E> rather than a double- 
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minimum landscape in Figure [2]A Knowing that the 
growth of the distance lj = — rj-±\ is a crucial fac- 
tor in the merging of energetic minima, we can suppose 
that lj is a subject to a similar interplay as dj, namely, 
the repulsive cores unable beads j — 1 and j + 1 to close 
up, while the stretching of lj is still possible. With the 
sufficiently high noise, this asymmetry could lead to the 
rise in (I) and the domination of the single-minimum po- 
tential topology. This transition seems inherent to the 
system and it is present regardless of the spatial corre- 
lations in noise. Nevertheless, once the single-minimum 
state prevails, the bead can explore the well, provided 
there is enough relative forcing, and for A ^ this forc- 
ing drops dramatically. In result, despite high noise am- 
plitude, the bead is trapped near the minimum, so the 
angles between modules cannot vary as much as in the 
non-correlated case. This leads to narrowing of 
around i/jj = 180°. 

X. SUMMARY 

Summarizing our research, the most salient conclusion 
one can draw is that spatial correlations in thermal noise 



have an overall inhibitory effect on the system. This man- 
ifests in the general attenuation of the chain geometry 
dynamics both in the time domain, where the polymer 
tends to preserve its current shape, and in the temper- 
ature domain, where the evolution of statistical chain 
properties is retarded. It is also in agreement with our 
previous findings that the presence of nonzero correla- 
tions reduced the ability of chain to transfer between dif- 
ferent conformations [13]. 



Such a behavior is not a surprise, as we can perceive 
the introduction of the spatial correlations into the ther- 
mal bath as a freezing of environment, and, in the limit 
A — >• +oo, this should also lead to the complete atten- 
uation of the system dynamics. In this context, it is 
not solely the temperature of thermal bath that influ- 
ences the system behavior, but also the structure of en- 
vironment. Our approach, which decouples the temper- 
ature from the environmental correlation length, allows 
for more plasticity than the explicit simulation of ther- 
mal bath particles, yet requires proper scaling to avoid 
unphysical situations. 
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